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This work describes the algorithm and encloses the program (language 
"Fortran-1" ) for identifying segmerits ‘of information obtained, from the 
telescopespectrometer "Filin", t^at are "suspicious" for the presence 
of an x-ray source. In .accordance with [1] the proposed algorithm is 
an algorithm of the lowest level. The information that is freed of 
uninformative segments. is evaluated by the .organizers of the experiment 
and subject to further,- processing with the involvement of algorithms 
of a ^higher level. 
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ALGORITHM AND PROGRAM FOR INFORMATION PROCESSING 
WITH THE "FILIN" APPARATUS 

by L. S. Guxin, V. S. Mokrov, Ye. I. Moskalenko ?nd K. A. Tsoy 

This work describes the algorithm and encloses the program (language "For- 
tran-4") for identifying segments of information obtained from the "Filin" tele- 
scope - spectrometer that are "suspicious" for the presence of an x-ray source. 

In accordance with [l] the proposed algorithm is an algorithm of a lower 
level. The information that is freed of uninformative segments is evaluated 
by the experiment organizers and is subjected to further processing with the 
involvement of algorithms of a higher level, 

1. Content of Experim ent 

On the station "Salyut-4" that was launched on 26 December 1974 the 

"Filin" telescope-spectrometer was installed that was designed to record 

radiation from space sources in the range 0.2-10 kev. In order to obtain 

spectral characteristics of the source the given energy range was divided into 

6 subranges: 0.2-0, 6 kev (channel Fl-5) , 0.6-0. 9 kev (Fl-6) , 0.9-2 kev (Fl-7) , 

2-3.1 kev (Fl-l), 3. 1-5. 9 kev (Fl-2), and 5*9-10 kev (Fl-3). In addition, there 

were two more summary channels: Fl-8 (0.2-2 kev) and Fl-4 (2-10 kev), and a 

channel that recorded the background of charged particles, Fl-12. The mechanical 

slit collimators limit the visual field of the x-ray detectors to 3x10° with 

respect to the half-width of the beam pattern. The instrument readings re- 

corded on magnetic t ape after primary processsing of information are values 
* 

Numbers in margin indicate pagination in foreign text. 



of the rates of counting in each channel. 


The "Filin" apparatus operated about 120 hours in different patterns. For 
a detailed description of the instrument and operating patterns see [ 2 ]} here vie 
will only examine the pattern of orientation in the local zenith at which the 
position of the optic axis of the telescope during movement of the station was 
maintained in the direction of the local vertical, while the visual field of 
the instrument was moved over the celestial sphere perpendicular to the broad 
side with a velocity of abort “ arg.min/s. Here the recording of signals from 
the observed sources is an .ojsceles triangle with width of the base about 9C s, 
which corresponds roughly tc JOO instrument readings in each of the examined 
channels. /4 

The abundance of obtained information significantly impairs isolation of 
weak sources on the background of charged particles. Therefore the task at the 
first stage is reduced to "rejection" of those segments of information where 
there are no sources. This work presents a method that makes it possible to 
isolate the places in the information in which a source is mathematically 
possible, and after verification of the physical criteria fed into the computer 
for the presence of a source to issue "suspicious" segments for more detailed 
analysis. 

The physical criteria include: 

1. The source in the examined operating pattern of the instrument is 
recorded in the form of a triangular -shaped impulse on recordings in different 
energy ranges, whereby the width of the triangle base cannot exceed 140 s, other- 
wise we have an anomaly, i.e., the passage through the metering devices of the 
instrument of a stream of charged particles. 

2. If the velocity of counting in the maximum impulse in the Fl-12 channel 
is comparable or greater than the corresponding amount in the channels that 
record the x-ray radiation, for example, in channel Fl-4, then "imitation" of the 
x-ray radiation by the charged particles occurs. 


2 


3. In sources of x-ray radiation the velocity of counting in the Fl-3 
channel must not exceed the value of counting velocity in the Fl-2 channel. 

This follows from the behavior of the curve of recording effectiveness of 
x-ray quanta [?], as well as ftom the fact that the sources that possess a 
sharply riiinr spectrum in this area until now were not observed. Therefore in 
such a case we will consider that we again have "imitation" of the source by 
the charged particles. 

2, Information and Algorithm 

The initial information before feeding into the program is preliminarily 
averaged for each channel such that the possible source is represented by 5-9 
points. In particular, in the examined operating pattern of the instrument the 
averaging was done for 48 values of the counting velocities. Segments of 
averaged information, where the instrument readings in more them 15 points in 
a row exceed 100 i/s are considered anomalous and are nullified by the program. 

Since we are interested in the segments which are "suspicious" if but for 
one of the channels that record x-ray radiation, then further discussions refer 
only to one channel. 

It is assumed that the result of measuring the random process is an addi- 
tive mixture of the signa] and the interference, and due to the discreteness of 
the measurements can be written in the following forms 

tfi m * Si ( 1 ) 


where Sj=f(t^) — signal and ^ — interference. 


Since it is known that the signal has a triangular shape, then can be 
presented in the forms 



0 

a(i- ) I* ‘ , 


( 2 ) 




3 



where 2k + 1— width of signal, a — its amplitude, and i 0 — time parameter. 


It is required that according to the measurements of Y^ the unknown 
parameters be determined (in our case a, i 0 , k) . With fixed values of i Q and 
k, the parameter a, using the least square method, is found from the expression 


2 


a ' pt*7 )&**'/) 2 (2« /3/ 


“here Jtf. (y f + V7M )/ 2 . 
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The segments of information for which a > C5^ v are considered "suspicious". 

Here C -Vfi . D — mean value of dispersion of measurements Y, with respect 
w av "^av av i 

to the session. The amount of D_„ is determined with respect to the background, 


*v 


whereby the background refers to the interval corresponding to 90$ of the area 
of the histogram constructed according to the measurements. In oefining D 

cp 

the following hypothesis was considered, following from the preliminary analysis 
of information: ^ H ' M — nu Nber of sources in a session, nj — width 


of source, N — total number of points) which means that the presence of sources 
does not noticeably distort the histogram of the background. 


The general algorithm looks as follows: 

1. With respect to the entire session for the values not exceeding 50 i/s 
we construct a histogram and define D 

av 

2. F'rom (3) for the assigned K at each point of the session we define a. 

3. The points in which a > 6" occurs and the criteria for the presence 

of a source as "suspicious" are fulfilled, are separated for more detailed analysis. 

_3_. P rogram and Results of Calculation 

The proposed algorithm was realized for the YeS computer in the 
"Fortran-4" language 
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The program operates with averaged initial information of the session 
divided into IN files SDSK(4lO). In the proposed text of the program, the / 7 

Information counting occurs from the magnetic disk ( YeS computer) . 


The SDSK structure t 

SDSK(l) — number of points in the file; SDSK(2) — number of session; SDSK(3- 
5) — date of session. Starting with SDSK(IO) the information referring to the 
given point I of the file is recorded in groups of 10 numbers each ( I< 40) i 
SDSK(lO+(l-l)xlO) — Moscow time, for example, 233409 — 23 hours, 34 minutes, 9 seconds; 
SDSK(ll-19Kl-l)xl0) --values of counting rates in channels Fl-1, . . . ,Fl-8, Fl-12. 

The groups of "suspicious" points between which the distance does not 
exceed K are stored in the file MAP. If this condition is violated, then the 
given segment together with the K-environs axe issued in print (subprogram 
PAIPER) . 

In the subprogram MID with respect to the segments of "empty" information 
the mean values are defined for the counting rates and the mean square deviations 
(j-1 3,12'. 

As a result of the work the programs are issued in print (for illustration 
we will limit ourselves to an examination of channels Fl-1, 2, 3, 12) « 


1. Initial information of session! 


time 

Fl-1 

Fl-2 

Fl-3 

Fl-12 

Number 
of ses 

1-31 

04.6 

203.9 

5*3.1 

*1.3 

98 

133* 

S3.* 

17V. 2 

487.* 

53.9 

99 

1321 

70.3 

1 4 * . C 

404.4 

33.4 

100 

139* 

*9.3 

143.1 

348.1 

38.3 

id 

12*1 

76.3 

150.2 

334.5 

* 7:1 

102 

1237 

SS.1 

15V. 3 

334.4 

79,0 

10» 

1222 

9S.*< 

174. 8 

321.4 

89.2 

10* 

1207 

1 09.* 

1*7.9 

318.7 

105.8 

103 

1U2 

1 09 • 3 

193.5 

321.8 

124.8 

10* 

1137 

120.3 

20>.4 

304*1 

130.1 

10? 

1122 

107.2 

1*7.5 

2*5.* 

115.4 

10® 

1107 

*7.7 

124.7 

1*2.6 

*0.2 

10’ 
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time 

Fl-1 

Fl-2 

Fl-3 

Fl-12 

Number of 
point of session 

loss 

44,4 

• l .6 

143*6 

17. 1 

110 

loss 

34,6 

* 0.1 

117.1 

51 .3 

111 

1023 

26,6 

62.7 

106,0 

50.4 

112 

1006 

26,1 

56.4 

6 $ . 7 

50.0 

1U 

953 

20,0 

44.2 

72.6 

43.2 


936 

14.2 

31.4 

S5.7 

39.6 

113 

923 

11,2 

24.3 

41.1 

36.2 

116 

90S * 

9.6 

23.6 

37.1 

36j,6 

117 


2. The histogram for determining the background, number of points by 
which the background is defined, and mean value of the background according to 
the session and 


i/sec 

Fl-1 

Fl-2 

Fl-3 

FI 

1 . 

*.0 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

0.0 

3 

o.o 

o.o 

0.0 

0.0 

4 

o.o 

0.0 

0.0 

0.0 

S 

1 .0 

0.0 

0.0 

0.0 

6 

o.o 

o.o 

1 .0 

. 0.0 

7 

42.0 

1.0 

1.0 

0.0 

8 

65.0 

3.0 

16.0 

0.0 

9 

13.0 

27.0 

39.0 

0.0 

10 

20.0 

52.0 

37.0 

0.0 

1 1 

9.0 

20.0 

12.0 

1.0 

12 

2.0 

5.0 

4.0 

0.0 

13 

2.0 

3.0 

6.0 

0.0 

14 

1 .0 

6.0 

1.0 

0.0 

15 

r.o 

7.0 

1.0 

1.0 

16 

0 .0 

1 .0 

2.0 

0.0 

17 

1.0 

4.0 

6.0 

7.0 

18 

o.o 

1 .0 

7.0 

21 .0 

19 

o.o 

0.0 

1.0 

31 .0 

20 

1.0 

1.0 

1.0 

32.0 

21 

0 .0 

3.0 

0.0 

7.0 

22 

o.o 

6.0 

0.0 

6.0 

23 

o . o 

7.0 

2.0 

4.0 

24 

1 .0 

4.0 

3.0 

0.0 

25 

0 • o 

1 .0 

2.0 

2.0 

26 

1 .0 

1 .0 

5.0 

0.0 

27 

O.o 

0.0 

2.0 

3.0 

28 

0*0 

1 .0 

0.0 

5.0 
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'se c 

Fl-1 

Fl-2 

Fl-3 

FI- 12 

29 

2*0 

0.0 

1.0 

1*0 

30 

0*0 

o.o 

0.0 

2.0 

31 

O.o 

1.0 

2.0 

3.0 

32 

o.o 

0.0 

1.0 

1.0 

33 

0 • 0 

0*0 

1 .0 

o.o 

34 

o.o 

1.0 

0.0 

1.0 

3 S 

1 .0 

0.0 

0 . 0 . 

2.0 

36 

o.o 

0.0 

1.0 

1 .0 

37 

0*0 

0.0 

1 .0 

3 .0 

36 

0.0 

0.0 

0.0 

Iff 

1.0 

39 

0.0 

0.0 

0.0 

9.0 

40 

0.0 

0.0 

0.0 

3.0 

41 

o.o 

0.0 

1.0 

0.0 

42 

o.o 

0 . 0 * 

0.0 

2.0 

43 

1 . 0 

0.0 

0.0 

3.0 

44 

1.0 

1.0 

1.0 

1.0 

4 $ 

o.o 

0.0 

0.0 

0.0 

46 

0. 0 

0.0 

0.0 

1 .0 

47 

0.0 

0.0 

0.0 

1.0 

46 

o.o 

0.0 

0.0 

2.0 

49 

0.0 

0.0 

1.0 

1.0 

90 

b.O 

0.0 

0.0 

2.0 


192 

127 

110 

11 0 


4.1 

40.0 

0.1 

12.2 


0.126 01 

o. 17E oi 0 

.126 01 

0.176 01 

i number 

of "empty" 

points between 

the "suspicious 

" segments, 

of counting rates 

in them and o . 

3 

0=1.2 8.12) 

• 

.ues of 

amount a in 

the "suspicious" points (for comparison o 

session 

is printed 

over them) and 

the segment of 

initial infi 


12 


av 


accor- 


matlon corresponding to these points together with K-environs: 


1 
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time 

Fl-1 

FI -2 

Fl-3 

FI -12 

number of j 


2 

2 

2 

9 



* 24.0 

244.2 

29*0 

27.2 



0 . 1*1 0 * 0 

. 11 * 02 

0.041 01 

0 . 12 * 02 



1*2 

1.7 

1.2 

1.7 


294744 

• 22*1 

214.0 

20.4 

42.0 

209 

2 S * 7 )J 

1211*0 

42?. 4 

70.0 

4 *.* 

204 


1010.4 

220.0 

47.0 

40.4 

20 * 


22*0 

24.4 

12.1 

21.1 

200 


141.0 

101,1 

22.7 

20.4 

201 


4 Q 4.7 

2 * 4.2 

12.1 

40.2 

202 


1117.0 

440.4 

• 2.4 

72.2 

209 


1470.1 

* 4 «.l 

03*7 

04.2 

204 


1244.2 

4 * 4 .* 

74.4 

7 *. 2 

202 


740.0 

* 0 ?.* 

4 J .2 

47.7 

20 * 


220*0 

221 .* 

12*2 

2**1 

207 


IS . 2 

20.2 

4.2 

1*.0 

20 * 


12 

12 

12 

12 



■ 0.0 

9.2 

40.0 

21.7 



0.741 01 0 

. 12 * 01 

0 . 14 * 01 

! * •. L4JL9-1- 


From 

an analysis of 

the trajectory data it 

follows that the 

"susp. clous' 


segment given as an example corresponds tc the source SCOX-1 that falls at the 
given moment in time in the visual field of the instrument--the brightest x-ray 
source in the celestial sphere. 


'KKJfN'AL PAGF. L 

"F POOR Gl’Aixn 


8 



Block Diagram of Program 


input and printing of initial information 
of session, nullification of anomalies, con- 
struction and printing of histogram, deter- 
mination according to it of mean background 
and dispersion according to session, and 
issuing of these amounts in print (subpro- 
gram "mid") 




Q 


to start cf cvcle 


Taxi of Program 


0001 


01**MIt>* tOSKUlo* .M10>» 

144<10> #?*<10> #*0<10> #*l(10> * 

214 <1 9 > #a«<io> .?A4U'09«10>. 
3TCH<20»10>#A2C200.11> #**? 0 0.41) 

0002 


OEM*! Ml* ?<?00#1**0’l#IM) 

000) 


421*10 

0004 


4Hl«? 

ooos 


442 H <3 * 1 0?) 

0000 


*AJU<1.104> 

000? 


NOT*© 

0.000 


00 10 2*2# 10 

0000 


Ml(i>*0 

0010 


00 10 J*1#41 

0011 

10 

HA?vI , J)«0 

0012 


00 *4 4*1 # 2 

0013 


00 11 1*1 # 10 

C0 14 


tQUXO. 

001S 


AC2)*0 , 

001* 


44 ( 1 ) *0 

001? 


00 11 J *1 , 30 

0010 


4A4 < J , t ) *0 • 

0010 

11 

com uui 

0020 


I4**21 

0021 


00 14 4*1,442 

0022 


ilAet? » 141 (SOS* <l) #2*1 #*10> 

0023 


2N*14«1 

0024 


4T«lMX(fet4(1M 

0021 


00 14 *1 *1 # 4T 

002* 


l ? l*“1 > *9 , *9# ?0 

002? 

*t 

M*T*4?T«1 

0020 


H*9*C*1*1?*10*1 

0020 


*T*1MX(S0S4C11>> 

0030 


23*11*1 

0031 


14*11*9 

0032 


041 11(3* 1|2)MTt <S0**<2> * 2*13 #2*1 « HOT 

0033 


00 *9? 2*2*10 

0034 


I2*t f *2-1 

0033 


23»«*»AMI ,l>»* 

003*r 


1 M*CIK( 22 >-100)491 #094 #494 

003? 

*01 

2M*9(t) -13)492 #493 #093 

0030 

002 

M?< 2 # 23) *0 

0030 


NS < 1 > *o 

00*0 


CO 10 49? 

0041 

093 

24*23*1 

0042 


4*<1>*0 

0043 


NA?<2,24)*4?T-1 
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0044 

0043 

0044 
0047 
004« 

0049 

0050 

0051 
0012 

0053 
0034 
0055 

0054 
0057 
0034 
0039 
0060 
0061 
0062 

0063 

0064 
0063 
0066 
0047 


NAP ( X * 1 >*PAPCX»1>*1 
60 1C 697 

694 !F(*$<X))69S«49S«494 
693 MAPI! «I3)»NPT 
494 DS(t)iNI(t)H 
497 CONlUUS 

70 CONlXkUC 

DO 14 I«2,10 
|f 

6*8048 (ID 
(6)14*14*12 

12 X9(*-1>71 .71 .73 

71 J»XUX<S»0 .3) 
19(4-30)72*72*14 

72 PAA(4*I)*9AA(JiXJ*1 
_60_10_14 

73 CONTINUE 
JP(3-RB(I>>13, 14,14 

13 KRW)«ICR(I)*1 
86 ( i > >$0 ( I ) *8*1 
A(X>*A(I>*4 

14 CONI XNUE 
19(6-1)74*74*43 

74 MAI I E (3 * 104) 


0046 

0069 

0070 

0071 

0072 

0073 

0074 

0073 

0074 
0077 
0074 

0079 

0080 
0081 
0082 

0083 

0084 

0083 

0084 
0087 
0086 • 

0089 

0090 

0091 


DO ( 5 J*1 , 50 

73 bRHE(3»1C3>.t,(PAR(J*IJfX**'10> 
. wftl I E (3 - 1 C*> 

DO •* 1*2*10 
J1*1 

oo n j *i 1 50 

X4(7ARCJf X)-PAil(Xi * I) >77.77*76 

74 I1»J 

77 CONI XNUE 

A<I)"MRCX1*X) 

X 2*3 0 

X9(A(X>>81 .41*74 
76 oO 60 |»I1 , SO 
S«P*RiJ'X)/A(f> 

19 (3-0. 03)79-79.60 
79 J2»4 

60 1C 81 
60 CONI IMJE 
41 R8U)»9L0A1 <X2> 

62 C0NIX4UE 
WRU E (3 * 109) 

CO 1C 84 

63 CHI PI0(KR*A,S6) 

84 CONI XNUE 
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original y-..cr- ... 
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0092 

0093 

0094 

0095 
009* 
009? 
009? 

0099 

0100 
0101 
0102 
OIOS 
010* 
0105 
010* 

' 010 ? 
0101 
0109 

one 

0111 

0112 

011S 

Oil* 

0115 

Oil* 

0117 

one 

0119 

0120 

8U5 

0123 

012 * 

0123 

.0126 

0127 

012 * 

0129 
0 1 3 C 

0131 

0132 

0133 

0134 

0135 

0136 
013? 
0131 
0139 


DO «5 X * 2 » 1 0 
83 RB(i)*Sfi(X) 

00 60 K a S » 9 • 2. 

*PT«0 

XN«*21 

K2«U«1)/| 

K1 a »2*1 

ICT a Q 

*1*0 

N 2*0 

*3*1 

N **0 

*5*0 

W 6*0 

00 13 1*2.10 
NS C I ) • 0 
SG<1)*0. 

15 SR<X>*0. 

DO >0 M*1,MC2 

REAOt?' XNMSDSK1X) » 1»1»*10> 
IN ■ X * ♦ 1 
00 19 1*1.6 

19 msu>*ifxx(sdsk<x>> 

NT*? $ < 1 > 

IF C7-D20.20.21 

20 ms<d«k 

WRXlE<3,1C6> 

WRHt<3,1C0) <HS(I),I»1,3) 
WRI IE <3# 101 > <H*U) » !■*#*> 

wrueo.io*) 

URllE(3.1C/> 

WRX I E {3 . 10*) 

21 CON I I RUE 

00 »0 HIsl.NT 

npt»npt*i 

J1"»-1 

00 *2 J«i I J1 

J2*J*1 

00 1 2 I«1 . 10 
22 TCHlj, I)»TCM(j2fI) 

25*104 (M-1)*1C 
TCHlk* 1)*S0SK<XS) 

00 23 1*2.10 
X1»X54l-1 

XFt74P(X > i))23*23«221 
221 X 3*7 XP ( X * i ) 

DO 225 J 1 ■ 1 # X3 
I4*24(J1-1)*2 





0140 

0141 

0142 

0143 

0144 

0149 
0144 /• 
014? 
014t 
0U9 

0190 

0191 
0152 
0193 

0154 

0155 
0154 
015? 
0158 
0199 
0140 
0161 
0142 

0163 

0164 

0165 

0166 
016? 
0165 

0169 

0170 

0171 

0172 

0173 

0174 

0175 
0174 
0177 
017* 
0179 

0150 
0161 
0182 

0183 

0184 

0185 
0184 

0187 

0188 
0189 


XFUP?-MP(X*X4>>23*222*223 
222 SDS»<X1)a-9. 

60 <0 23 
229 X 4* 1 44 1 

XFUPT-PAPU ,X4) >222*222. 225 
225 C0NM5UE 

23 TCHUC.I>*SDSK(X1> 

5T*5T*1 

XFCM-K)Sfi.24.24 

24 5f*5T*1 

00 24 X*2 . 10 
KftUUO 
DO 26 J"1 .8 
XPOCH(JiX>>25*25*26 

25 KR<l>»K*(X>*1 

26 COMMPUE 

00 51 1*2.10 
XF<RR(X>>26. 28*27 

27 A(X>«*0.001 
60 10 31 

28 A<X>*0. 

F0N*<TCH(1 . X ) *TCH <8 * X) ) 72 

DO 29 J *51 .5 

J1*5*i-J 

» I)4TC«(J1 . X)-2*f04> 

29 A<X>*A<1>*3»S/<51«8> 

X F 1A ( I > >36 « 301 *301 

30 A < X > «0 • •’ 

60 1C 31 

301 J 4*6 

00 504 J3*1 .52 

XFUCHtJ3.I>-TCH(81 * X> >304*303*303 

303 J4*24a1 

304 -CONI I5UE 

1F(J*.-52>305*30.30 
309 J 5*0 

53*81*1 

DO 507 J3*53 • 5 

IFUCN(J3.1>-TCH(51»X)>307*306*306 

306 j 5* J 5* 1 

307 CONI I MJE 
IF(25-K2)31 .30*30 

31 CONIXNUE 
A<1>*TCHU1,1) 

J ■" 1 

00 513 * 1*2*10 
XF(8(X>-R81X>>313*311.312 

311 J-0 

60 10 313 

312 2*1 


| 


5 


13 




0190 

313 

01«1 


0192 

32 

0193 


0196 


0195 


0196. 

321 

019? 

322 

0198 


0199 


0200 

323 

0201 

324 

0202 

32$ 

0203 


0204 


020$ 

32* 

020* 

33 

020? 

36 

0208 

3$ 

0209 


0210 


0211 


0212 


0213 


0214 

3* 

0215 


021* 


0217 

37 

0218 


0219 


0220 


0221 

371 

0222 


0223 


0224 

372 

0225 


022* 


022? 


0228 

373 

0229 


0230 


0231 

374 

0232 


0233 


0236 


0235 

38 

023* 


0237 

381 

0236 

382 


C( to f Hwf 

IF(*!>33. 33 *32 
S-A 1 2 ) 

00 922 1-2/9 
V»A1.' j 

XF(V»S)322«322*321 

S*V 

CONI HUE 
$ 1 * A (10) 

IM»-S1>323,323*32* 

J*0 

X F (<f > 33 * 33 *325 

S»MJ) 

81*AC4> 

XF(6-I1>32*i32**33 

J«0 

XP(j)451 *361 *3? 
XF(61)384. 384*35 
CALL PX0(6S*SR* 86) 
WR1 1 E (3 * 104). 

ni«6 

DO 36 1-2,10 
NS(i>*0 
SR ( 1 ) -0 • 

SOCD-O. 

R4*0 

60 1C 386 
j.NFT-ICHI 

N**6*41 

AZ(6*,11)-FL0AT(J) 
DO 3 ? 1 1-1,10 
AZ (6* , X ) -A (X ) 

R3«-A2 

13(66)372*372,376 
DO 373 J1*1*K2 
X1*6PT»IHJ1 
PAR V J 1 ,1)-FL0AT < X 1 3 
00 373 1*2*10 
PAR(J1*X>«1CM(J1*X> 
R4*1 
N2**2 
62*62*1 
J-NM-K2 

PAR(K2*1 {-FLOAT <J) 
DO 38 1*2*10 
PARU2* I)-1CH(K1 ,X> 
60 1C SO 

17(63)382*383*346 
N3 *6 3* 1 


I 

r 

t 

* 


I 


i 

0 2 51 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 
0286 
0281 
0282 
0283 


0284 

0285 

0286 
0287 


0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 


60 10 374 

383 1F(46}452i4S2.34 

584 W*I1E<3#103><R8<I>#I»2. 10) 

CAU PAIPERU2«PAR#42#K6> 

WR1 1 E<3 r 104) 

389 N2»« 

R6»0 

R4*0 

H5**2 
60 10 50 

386 j F <45)41 «41 #387 

387 N5»A5-1 

60 10 50 
41 N1»1 

DO 45 1*2(10 
1F<1CH<K.I>)6^»45#44 
44 NS<l)*NS<I)*i 

S6<D"S6<I>*TCH<Kf2>*-KH(K«2) 

sro>*smi>*tcmu#i> 

4» CON1INUE 
60 1C 50 

451 IF <40)366,386, 381 

452 q0 456 11*1#N2 
DO 454 1*2 » 1 0 
TF<*’AMI1,i))454#454#4S3 

453 MS < 1 ) <NS ( X ) +1 
SR<1)*SR<!)*PAR<I1 »X) 

SC U) <S6<D4PAR<X1 #D»PAR<I1 , J) 

454 CON1X4UE 
N1 *1 

60 1C 395 

50 CONI X AUE 
IF<M)S2, 32,51 

51 CAU p ID (AS , SR * SP) 

GO 10 60 

52 CALL FAIP£H<A2,PARf N2,N6) 

60 C0N1IAUE 

100 FORP AT (3X5) 

101 FORPAK3X3) 

102 F0RPAT(X7,1X,9F10.1 #17) 

103 FORPAT<8X,9F10.1> 

104 F0RPAT<1H ) 

105 FORPAT<I7,1X,9Fio.1> 

106 FORPAK////) 

107 FolPAT <1X, ’ TIME' *SX. * Fl-1 ' ,6X« ' FL«2» ,6X, * FL-3* 
1 ,6*» ' FL-4' ,6X, ' FL-5' , 6 *# ' Fl-6‘ #6X. ' Fl»7' 
2,6X,*FL*8' ,5X« * F L»1 2 • »4X, r NPT ' ) 

10* F0RPAT</' histogram to define background '#/) 

109 F0RPAT</» *OHI' i »H0.1./) 

0E8L6 SUBCHK 
ENO 
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iiisMiiii: 


0001 

0002 

0003 

000* 

000 S 

ooo* 

0007 

0001 

0009 

0010 
0011 
0012 
0013 
001 * 
0013 


SUBHClTINE PAXPER<AZiPAR,N2iN*> 
DIMENSION PAR(200*10>f AZCZOOfll) 

00 1 J1*1 .*6 

. MT«1F IX(AZIJ1 f 1 > > 

J*IP XX (AZ (til ,11>> 

1 UAXIE(3»102>MTf <AZ(Z1f X)«X a 2#lO)iJ 
WRX I E (3 * 1 0*) 

DO * Z 1 *1 « 1*2 
J*XMX(PAR(X1,1>> 

Z WR1 1 E (3 > 1 0*) (pAR( XI •!) »X«2» 10>*J 
10* FORPAT(aX'9Flo.1'X?> 

10* FORPAKIH ) 

102 F0R*AT<X7i1X#9H0.1*X?> 
return 
END 


M 


C Co*. 


SUMoUriNC rtj • «*R ♦ 36) 

0002 

• 

DIHem$ION W!j(1) ,SR(1) *80(1) 

0003 


DO 3 1*2*10 

000* . 


XF(*is<I)-1>*r'>l2 

0005 

* 

SR<X)»-0,01 

0 Oc 6 

1 

SG< X > "0* 

0007 


CO To 3 

0006 

2 

S*^S<I)*SC(d-SR(X>*S*(X) 

0009 


S*S/<NS<I>*<11S<X*-1>> 

0010 


SG<X>*SQRT(S» 

001 1 


SR(n*SR(X3/DS(il 

0012 

3 

CONTINUE 

0013 


WRXTE<3'105>(N5<|) U 8 2M0> 

0014 


WRXtu< 3»103> (SR<f> U*2*1C> 

OOi $ 


WRXTe<3'107><SG<I>M«2MO> 

0016 

103 

F0RMaT(6Xi9F10»iX 

0017 

107 

F0RnaT(6X«9E10|2) 

0016 

103 

FORNaT(6X»9X10) 

001 9 


RETURN 

0020 


END 



v * "t - 

K ‘ 
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